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Abstract 

We derive the virial theorem of the relativistic two-body system for the study of 
the S-meson physics. It is also shown that the solution of the variational equation 
always satisfies the virial theorem. From the virial theorem we also obtained 
fil = — Ai = (p 2 ) = 0.40 ~ 0.58 GeV 2 , which is consistent with the result of the 
QCD sum rule calculations of Ball et al. 
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I. Introduction 

The U-physics provides many important and attractive physics informations, 
and is under active experimental and theoretical investigations. The 5-factories 
at KEK and SLAC will give valuable clues for the better understanding of the 
Standard Model (SM) and the beyond. We expect the 5-meson system will show 
the CP- violation phenomena |l|], for which we have only the — > tttt decay 
for more than 30 years. The mechanism of CP-violation through the complex 
phase of the Kobayashi-Maskawa three family mixing matrix || in the Weinberg- 
Salam model is presently the CP-violation within the SM. In order to understand 
and precisely test the SM, it is essential to know the values of the Kobayashi- 
Maskawa matrix elements, especially through the determinations of V u h and 
Vtd, and confirm the unitarity triangle both of which will be best probed at 
the forthcoming .B-factories. 

Recently, it has been an important subject to obtain an accurate value of 
the kinetic energy, /i 2 (= — A x = (p 2 )), of the heavy quark inside 5-meson in 
connection with the heavy quark effective theory (HQET) || ||. Ball et al. @ 
calculated ^ 2 using the QCD sum rule approach and obtained /i 2 = 0.50 ± 0.10 
GeV 2 for S-meson, while Neubert § obtained -Ai = 0.10 ± 0.05 GeV 2 . Neubert 
also derived || the field-theory version of the virial theorem within the HQET 
framework. Based on the theorem, he implied that the result /z 2 ~ 0.5 GeV 2 of the 
QCD sum rule calculations of Ball et al. is too large. However, it should be noted 
that Refs. and || differ in the choice of the 3-point correlation functions used 
to estimate the matrix elements of interest. The difference in the numerical values 
obtained in these two calculations is understood in terms of the contributions 
of excited states, which in principle must be subtracted in any QCD sum rule 
analysis. In practice, this subtraction can only be done approximately. Therefore, 
the numerical differences between Refs. and || indicate the limited accuracy 
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of the QCD sum rule approach. Bigi et al. [10 derived an inequality between 
the expectation value of the kinetic energy of the heavy quark inside the hadron 
and that of the chromomagnetic operator, (p 2 ) > |(My 2 — Mp 2 ), which gives 
/i 2 > 0.36 GeV 2 for B- meson system. However, Kapustin et al. ]TTJ showed later 
that this lower bound could be significantly weakened by higher order perturbative 
corrections. Previously we also calculated |12[ the value of (p 2 ) by applying the 
variational method to the relativistic Hamiltonian, and obtained (p 2 ) = 0.44 GeV 2 . 
Similarly de Fazio [[U| computed the matrix elements of the kinetic energy operator 
by means of a QCD relativistic potential model, and found /i 2 = 0.46 GeV 2 . 

Besides the above theoretical calculations of /^ 2 , Gremm et al. [14] extracted 
the average kinetic energy by comparing the prediction of the HQET |L5[ with the 



shape of the inclusive B — > Xlv lepton energy spectrum [|nj for Ei > 1.5 GeV, in 
order to avoid the contamination from the secondary leptons of cascade decays of 
b — > c — > slu. They obtained — Ai = 0.19±0.10 GeV 2 . Combining the experimental 
data on the inclusive decays of D — > Xeu, B — > Xev and B —>■ Xtu, Ligeti et 
al. JT7| derived the bound of /i 2 < 0.63 GeV 2 if A > 0.240 GeV, or /i 2 < 0.10 
GeV 2 if A > 0.500 MeV. Li et al. |1| obtained the value of -X x centered at 0.71 
GeV 2 from the analysis of the inclusive radiative decay B —>■ Xg'-f [19] within the 
perturbative QCD framework. Related with the comparison of various theoretical 
calculations of /i 2 , we note that Ref. [[2(J emphasizes that one has to be careful 
when comparing the values of — Ai obtained using different theoretical methods. 
For instance, QCD sum rule determinations of —X\ for the ground state heavy 
mesons and baryons are affected by a renormalon ambiguity problem |2(J . 

Concerned with the phenomenological importance of the numerical value for 
the kinetic energy of the heavy quark, we would like to derive the virial theorem of 
the two-body system within the relativistic potential model approach for the study 
of the .B-meson system. We show at the end that /z 2 ~ 0.50 GeV 2 of Ball et al. JF| 
is consistent with our virial theorem result. In Section II we present the variational 



3 



analysis of the relativistic quark model for 5-meson system, and compare it with 
the HQET. In Section III we derive the virial theorem of the relativistic two-body 
system which is appropriate for 5-meson system, and we show that the solutions 
of the variational equation always satisfy the virial theorem automatically Section 
IV contains discussions and conclusions. 



II. Relativistic Quark Model 



For the study of the bound state properties of hadrons which contain both 
heavy and light quarks like 5-meson, it is appropriate to use the relativistic quark 
model |12 , plf| , which is a potential model approach with relativistic kinematics. 



In the relativistic quark model, the Hamiltonian for 5-meson is given by 

2 

H = \Jv> 2 + M 2 + yp 2 + m 2 + V(r) » M + ^ + yp 2 + m 2 + V(r) (1) 

in the 5-meson rest frame, where M and m are heavy and light quark mass 
respectively. For the potential in (|1|), we use the Cornell potential |2l|, |22|, |23|, 

V(r) = + Kr + Vo, a c = -a s , V = constant . (2) 

r 3 

It is difficult to solve the eigenvalue equation of the Hamiltonian operator ([!]), so 
we use the variational method. The variational method is particularly useful in 
the present work, since we will show in the following section that the solutions of 
the variational equation always satisfy the virial theorem. 

In the variational method, the expectation value of the Hamiltonian is calcu- 
lated with some trial wave function which has a variational parameter. The value 
of the variational parameter is determined by the stationary condition (variational 
equation or so-called gap equation). We take the variational parameter which has 
the dimension of mass. Then from the dimensional analysis, the expectation value 
of each term in the Hamiltonian can be expressed as 

(p 2 ) = CV, 



(\/p 2 + m2 ) = + a 2 m 2 //i + 0((m 2 //i) 2 ), 



a 



(-y + Kr + Vo) = a c (-hii) + K(b 2 /fx) + V , (3) 

where \i is the variational parameter of mass dimension, and C, czi, a 2 , &i arid 62 
are dimensionless numerical constants. Collecting the terms in (|||), we have 

C 

(H) lx = E(n) = M +[V + {a 1 -b 1 a c ) f i+{a 2 m 2 + b 2 K)/fi} + —fi 2 

= M+[\/ + ^ + 7 //i] + ^ 2 , (4) 

where 

(3 = d\ — b\a c and 7 = a 2 m 2 + 62^- (5) 
We neglected 0((m 2 /{i) 2 ^ terms in (f|). Then the variational equation reads 

^E{n) =(3- 7//X 2 + = 0. (6) 

Rather than solving this equation numerically, we obtain the solution as a power 
series in 1/M, since M is large. Noting that the solution is given by p, ~ for 
very large value of M, we expand 

fi = h + h 1 ^ + h 2 j^ + ---. (7) 



Substituting (|7|) into (0) and matching order by order, we get 

Then we get E(ji), the 5-meson mass Mb, as a power series in 1/M, 

E{fi) = M + (v + 2^) + °- + O(^). (9) 



M, 



Therefore we have performed the 1/M expansion for Mb in the relativistic quark 
model in the framework of the variational method. Then, let us compare the series 
expansion in (^) with that of the HQET |5|, |5j which is written as 

M B = M + A + JL(T + i/ B n) + O(^), (10) 



where vb = 1/4 and —3/4 for vector and pseudoscalar 5-meson respectively. The 
zeroth order term A in ( |T0| ) which is the contribution from the light degrees of 
freedom corresponds to the sum of the expectation value of the light quark kinetic 
energy and that of the potential energy, 



A < — > (^/p 2 + w 2 + V(r)) = V + 2^(3. (11) 

The heavy quark kinetic energy term T in fllOD has the correspondence 

T^{p*) = C 1 -. (12) 

If we had included a spin dependent potential which is inversely proportional to 
the heavy quark mass, the expectation value of which would have corresponded to 
the chromomagnetic interaction term in fllOD as 

Y^L^iy\ = Ll v \ with v s = — Sl -s 2 V 2 (-—), (13) 

where v s is the nonrelativistic spin-spin interaction potential. 

III. Virial theorem 

In classical mechanics, the virial theorem for a system with kinetic energy K(p) 
and potential energy V(r) is based on the relation 

|(r • p) = v • p + r • ±p = (^(P)) ' P + r • (-W(r)), (14) 

where K(p) is the kinetic energy given by K(p) = J F • dr = J v • dp. The time 
average of the left hand side of (|T4]) vanishes for a periodic motion or a bounded 



motion in an infinite time interval, hence we get the virial theorem, 

d d 
(p ■ -K-K(p) - r • — V(r)) thae average = 0. (15) 
op or 

This theorem holds for both nonrelativistic and relativistic kinematics which have 
the following relations respectively, 

p 2 OK p 2 

p = mv , K = g-, p. — = ?- = 2K- (16) 
2m Op m 
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mv r— dK p 2 

P = n i n K=yl-p l + m i , p- 



Vl - f 2 ' <9p \/p 2 + m 2 ' 

In quantum mechanics, the virial theorem for a system with Hamiltonian H = 
K(p) + V(r) is based on the relation 

-< r .p> = _<[ r .p,ff]) = _ (r .(-_) + (_).p). (17) 

For stationary states represented by eigenfunctions of H, the left hand side of ([1 
is zero. Hence we get the virial theorem written as 

dV dK 

<"-<-ar> + <ap->-p>-°- < 1! 



However, the functions which satisfy the virial theorem (|T8|) are not restricted 
to the eigenfunctions of H. The solutions of the variational equation also satisfy 
the theorem ([18]). With the variational parameter /i of the mass dimension, the 
expectation value of the Hamiltonian is expressed as 

E(fi) = (iP(fi)\H\^)) = (K(p) + V{r))», (19) 

for central potential V = V(r). In general, the kinetic and potential energy func- 
tions in (|T9|) can be expanded in the Laurent series as 



K(p) = £ k n p n , V(r) = £ v n r n . (20) 

n n 

The expectation values of each terms are written in terms of fi, 

(p n ), = a n fi n , (r n ), = b nf i- n , (21) 

where a n , b n are numerical constants which have no dimensions. The expressions 
in (pi] ) satisfy the relations, 

= {p t/ ] '" ^ ) - = - {t t/" ) - (22) 



From (|19h-(|22T) , we get 



"^ ( " ) = <? ^ w -^ nr)) - (23) 



The equation ( p3[ ) means that the variational equation, dE(/j,)/djj, = 0, is equiv- 



alent to the virial theorem ([18]). If we have used the parameter /z having the 
dimension of length, we would have gotten the same result except for the overall 
sign, which is irrelevant. 

We can also show the relations more explicitly. For the expectation value of 
the potential, 



(V(r))„ = / dTr | V(r; fi)\ V(r) =n I dr r 2 A(jir)V (r) , (24) 
where we defined A(fir) such that it satisfies 

j d 3 r | %Ij{y- /i)| 2 = J drr 2 \ R(r; /i)| 2 = // 3 J drr 2 A(fir) = 1. (25) 
Changing the integration variable in ( p4|) by x = /ir, 

(V(r))„ = J dxx 2 A(x)V(x//j) } (26) 

then we get the relation 

^(V(r))^ = fdxx 2 A(x){--V'(x/fi)} 

3 [ J_2 



/ drr 2 A(/ir){-ry'(r)} 
-<r^y(r)>, (27) 



For the expectation value of the kinetic energy, we can follow a similar procedure in 



the momentum space to get the desired result jj,-^-(K{p)) ^ = (pj-K{p)). Therefore 



we obtain the relation (53) again. 

In the above, we have shown that ip(fi) with Jx being the solution of the vari- 
ational equation, dE(fi)/dfi = 0, automatically satisfies the virial theorem given 
by 

(i>(jl)\ P —K(p)\i>(jl)) = (^)\ r -V(r)\^)), (28) 

even though it is not an eigenfunction of H. We emphasize that the solution set of 
the virial theorem contains the functions determined by the variational method as 
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well as the eigenfunctions of the Hamiltonian. Repko et al. []24| previously proved 
that variational principle implies the virial theorem, and Lucha et al. |25] derived 
the relativistic form of the virial theorem. Let us show explicitly that the virial 
theorem is satisfied for the system given by the Hamiltonian of ([!]) and (Q). 
For 5-meson which has the Hamiltonian ([!]) with K(p) = p 2 / (2M) + y/p 2 + m 2 , 
the virial theorem (E28|) is written as 



<S + V5fe > = (r ^ (r)> - (29) 

The left hand side of (p9|) , which is related with the kinetic energy, can be calculated 
from (E), 



'M ^ Vp 2 + ™ 2 Af VVP Vp 2 + ™ 2 



— — + (ai/i + a 2 m 2 /ii) - 2m 2 (aifi + a 2 m 2 / » 

M o{m z j 

M 



+ (a ly u — a 2 m ///). (30) 



The right hand side of (p9[), which is related to the potential energy, is also given 
from (||), 

3 n 

(r—V(r)) = (-£ + Kr) = a c (M + K(6 2 //i). (31) 
ar r 

Putting ( p0|) and (|3~ID in (|29|) and dividing by fi, we get the variational equation 
(0). Therefore we have shown that the virial theorem is satisfied for the solution 
of the variational equation of the Hamiltonian given by ([[]) and (Q). 



IV. Discussions and Conclusions 



We have shown that the relativistic quark model combined with the variational 
method can give many useful results. The solutions of the variational equation 
always satisfy the virial theorem automatically, and the 5-meson mass Mb is ex- 
pressible as a power series in 1/M. Combining fl29f) and fl3l|) , we get the relativistic 
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virial theorem within the Cornell potential, 

. f[ 2 ) = (r-W) 
M VP + m 



QJ 

(_£ + Kr ) = « c ( &l/i ) + K{b 2 /n). (32) 



With the usual input values |T], ||, [23| of a s = 0.21 ~ 0.36, K = 0.19 GeV 2 , m = 
0.15 GeV, we obtain /x 2 = 0.40 ~ 0.58 GeV 2 for Gaussian trial functions, which 
is consistent with the result /i 2 ~ 0.50 GeV 2 of the QCD sum rule calculations of 
Ball et al. J7|. If we use exponential trial functions, we get higher values of /z 2 . If 
we rather had used the virial theorem of the nonrelativistic two-body system 

<^ + ^> = <r-W> = <^ + ^>, (33) 

we would have obtained the value /i 2 = 0.10 ~ 0.15 GeV 2 , which is very small 
compared with the above relativistic result. 

If we apply the virial theorem to the system of one-body in an external poten- 
tial, it reads 

<^> = (^V(r)), (34) 
with the nonrelativistic kinematics for the heavy quark. This corresponds to a 
one-body (heavy quark) motion in a fixed external potential of the background 
system which produces the potential in the meson rest frame. Compared with 
the hydrogen atom, where the nonrelativistic one-body virial theorem (through 
the reduced mass) is also applied, the roles of the heavy (proton) and light (elec- 
tron) degrees of freedom are reversed. Within the potential model approach, the 
background system corresponds to a valence light quark and virtual gluons, the 
former carries the compensating momentum against the heavy quark motion and 
the latter contributes as a potential energy. Considering these, the correct virial 
theorem for 5-meson system within the potential model approach should be the 
form of the two-body closed system as given in (j32|). 

In conclusion, we derived the virial theorem for the relativistic two-body sys- 
tem, and numerically obtained /i 2 = — Ai = (p 2 ) = 0.40 ~ 0.58 GeV 2 , which is 
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consistent with the result of the QCD sum rule calculations of Ball et al. It is also 
shown that the solution of the variational equation always satisfies the virial the- 
orem. As our final comment, we note the recent observation |2D) on the mixing of 
the operator for the heavy-quark kinetic energy with the identity operator, which 
implies that the parameter Ai of the heavy-quark effective theory is not directly 
a physical quantity, but requires a non-perturbative subtraction. Concerned with 
the phenomenological importance of the numerical value for the kinetic energy of 
the heavy quark, it would be the most urgent to try all the possible attempts on 
the physical understanding of the kinetic energy of the heavy quark inside .B-meson 
system. 
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